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We study the effect of photothermal fluctuations on squeezed states of light through the photo- 
refractive effect and thermal expansion in a degenerate optical parametric amplifier (OPA) . We also 
discuss the effect of the photothermal noise in various cases and how to minimize its undesirable 
consequences. We find that the photothermal noise in the OPA introduces a significant amount of 
noise on phase squeezed beams, making them less than ideal for low frequency applications such as 
gravitational wave (GW) interferometers, whereas amplitude squeezed beams are relatively immune 
to the photothermal noise and may represent the best choice for application in GW interferometers. 



I. INTRODUCTION 

Optical squeezed states are used in many areas of 
quantum optics to improve the sensitivity of measure- 
ments to beyond the shot noise limit (SNL). For example, 
squeezed states can be used in interferometers Q , spatial 
and spectroscopic measurements Q, and potentially 
to improve the quantum noise limit of gravitational wave 
(GW) interferometers 0,0,0- Optical parametric am- 
plifiers (OPA) are often the systems of choice to produce 
squeezed states since, in theory, they can produce states 
with very high levels of squeezing. The level of squeez- 
ing that can be produced in these systems is limited by 
the introduction of noise from a variety of sources. The 
noise sources that have been reported to limit squeez- 
ing in these systems are pump noise H, 0] an d seed 
noise [13 [HI. 

In experiments reported to date, the maximum amount 
of squeezing inferred before detection is around 7dB ^(| • 
This result, and in most results, the maximum squeez- 
ing is measured at sideband frequencies above 1MHz, 
rather than at lower frequencies where ideally, greater 
squeezing is predicted. There has been recent interest 
in producing squeezed states at lower frequencies, pri- 
marily for use in gravitational wave detectors. For such 
states the squeezing bandwidth should cover the GW de- 
tection band (10Hz to 10kHz ). S everal results have been 
published below 300kHz 0ll![l4| with the lowest re- 
sult 280Hz These frequencies represent a different 
regime experimentally to the majority of squeezed state 
production and as such other potential limiting low fre- 
quency noise sources need to be considered. 

One such effect which is large at low frequencies in op- 
tical cavity systems is the photothermal-effect-induced 
noise 0, El- This effect has not previously been 
reported in an OPA system to our knowledge. By na- 
ture, photothermal effects are important at low frequen- 
cies and may be significant for limiting the generation 
of squeezed light in the GW detection band. This effect 
is investigated theoretically in an OPA cavity system in 



this article. 

The photothermal effect can be described as the ab- 
sorption of optical power in a medium causing a tempera- 
ture change to the medium. This effect may be significant 
in most nonlinear crystals since many have relatively high 
absorption rates. For example, the crystal MgO:LiNb0 3 
has the linear absorption rate of 4%cm~ 1 at 532nm [l8| . 
High absorption rates coupled with the high circulating 
power required for strong nonlinearity result in a large 
amount of optical power absorbed into the crystal, which 
may cause a significant temperature change. The aver- 
age temperature change due to the power absorbed in 
the crystal can be compensated for by using a temper- 
ature controller, and does not pose significant problems 
for most experimental systems. 

Instead, we focus our investigation on the effect of pho- 
tothcrmally induced temperature fluctuations, caused by 
fluctuations in the circulating power in the OPA cav- 
ity. The circulating power fluctuations could have classi- 
cal and quantum origins, or in the case of a shot-noise- 
limited system, only quantum mechanical origin. The 
photothermal noise caused by thermal-expansive noise 
and thermal-refractive noise has two degrading effects on 
the production of squeezed light in the OPA. The first 
effect is via fluctuations in the nonlinearity. This arises 
as the nonlinear strength is temperature dependent due 
to the phase matching condition. The second effect is 
via optical path length fluctuations. The temperature 
fluctuations cause the optical path length to change, po- 
tentially causing a detuning of the optical cavity. These 
effects appear as l/(^ 2 + &r) in variance, where fi-r is 
the thermal relaxation cutoff frequency of the nonlinear 
medium, and are therefore primarily significant at low 
frequencies. 

This article is organized as follows: In Section [HJ we 
write down the equations of motion for the fundamental 
and second-harmonic modes in an OPA with extra terms 
required to take into account the photothermal effect. In 
Section IlIII we quantify the fluctuating photothermal ef- 
fect. In Section Till Al the relation between fluctuations 




FIG. 1: A Schematic of the OPA cavity. Ai n and Bi„ are the 
input fields to the OPA cavity, A out and B ou t are the output 
fields, and a and b are the intra-cavity fields at the fundamen- 
tal and second-harmonic frequencies respectively. 7a™b, "/abi 
7^ c 6 , and 7^7/ are the cavity damping constants associated 
with the input- coupling, output-coupling, intra-cavity scat- 
tering, and intra-cavity absorption respectively. The index 
(a,b) refers to the fundamental and second-harmonic frequen- 
cies respectively. Sv°^, 5v°^ b , and Sv^ £ are the associated 
vacuum fields that couple in. 



in the power absorbed into the crystal and in its temper- 
ature is described. In Sections IIII Bl and IIII CI the cou- 
pling of the temperature fluctuations to the fluctuations 
in the nonlinear coupling strength and cavity resonance 
frequencies through the photo-refractive effect and ther- 
mal expansion of the crystal is described. In Section ll VI 
the equations of motion with these photothermal contri- 
butions are solved, and the quadrature field amplitudes 
in both amplitude and phase quadratures are studied. In 
Section [V] the squeezed/anti-squeezed quadrature vari- 
ances with the inclusion of the photothermal noise are 
derived and plotted. In Section IV Al we discuss the re- 
sults for standard experimental parameters and for the 
shot-noise-limited case. In Section IV Bl we consider the 
effect of squeezing with the photothermal noise on grav- 
itational wave interferometers at low frequencies. The 
conclusions of the paper are summarized in Section lv"fl 



II. FIELD EVOLUTION EQUATIONS IN THE 
DEGENERATE OPTICAL PARAMETRIC 
AMPLIFIER 

In this section, the model of the OPA based on the 
Heisenberg equations of motion is introduced, then these 
equations are linearized and additional terms for the pho- 
tothermal fluctuations are introduced. This sets up the 
formalism to include the photothermal fluctuations which 
are described in terms of the input fields in the follow- 
ing section. The modes can then be out-coupled and the 



variances calculated. 

Starting from the quantum Langcvin equation, the 
equations of motion for the intra-cavity fields at the fun- 
damental frequency, a, and at the second-harmonic fre- 
quency, b, are given by 

a = - K + iT ) « + z*^b + y^CA m e-^ 

+ V / 2^C* + K c + y/^F s 6v a a bs , (1) 

b = - {iujI + 7^) b - X -ea 2 + yfi&Bine-*-** 

+yfitf*8v? a + yffi? Sv s b c + sJ*t s K bs . (2) 

The fields and coupling rates here are shown schemat- 
ically in Figure ^ Ai n and Bi n are the fundamental 
and second-harmonic input fields to the cavity at fre- 
quencies u> a and Ub respectively (u>b — 2u a ). e is the 
nonlinear coupling constant. io c a and u b are the cavity 
resonance frequencies of the fundamental and second- 
harmonic fields. 7™ h , 7°1,*, 7^ c b , and 7° b h s are the cavity 
damping constants associated with the input-coupling, 
output-coupling, intra-cavity scattering, and intra-cavity 
absorption at both frequencies. dv^, and 5v^ b § 

are the associated vacuum fields that couple in. The fol- 
lowing commutation relations are satisfied: 

M = o, M] = i, (3) 

for s = A m ,B in ,8v° a ut 1 8v b Mt 1 8v s a c , 8v s b c , Svf s , and Svf s , 
and all others vanish. 

Transforming to the rotating frame of each field with 
a — > e luJat a, b — > e luJbt b, and similarly for the input fields, 
the equations of motion become, 

a = -( l Lu d a et + 1 t a ° t )a + e*^b+^2^A ln 

+ ^rt 5v OUt + ^C §< C + yf^Svf, (4) 

b = - (iuf* + 7^) b - l -ea 2 + ^CB ln 

+,fitf*5v? a + V^W K C + filF s 5v a b bs , (5) 

where the cavity detunings et = u c a — ui a and u> b et = 
u>l - LU b . 

The most common method of generating the analytic 
form of squeezed quadrature variances is to expand the 
operators about their steady state values and then lin- 
earize the resulting expressions to first order in the fluc- 
tuation terms |20|. To linearize the equations of motion, 
the following substitution for the annihilation and cre- 
ation operators are used: 

a = a + Sa, a' = a* + 5a' , (6) 
b = b + Sb, = b* + Sb\ (7) 

and similarly for the input fields A. in , B in . x for x = a, b is 
the complex expectation value (a;) and Sx is the operator 
for the fluctuations of x so that (Sx) = 0. 



2 



To obtain the equations of motion for the fluctuating 
components, we consider the following fluctuations due 
to the photothermal effect in the crystal: (1) fluctuations 
in the resonance frequencies of the cavity due to fluctua- 
tions in the length of the cavity (assuming that the laser 
frequency is perfectly stable), and (2) fluctuations in the 
nonlinear coupling constant due to fluctuations in the 
temperature of the crystal and in the length of the non- 
linear region. We then make the following substitutions: 



det - det , c det det - det 

(u)% et and co b et are real), 
e = e + 5e, e* = e* + 5e ] . 



- Su> 



del. 



(8) 
(9) 



The fluctuation terms are obtained by taking the fluctu- 
ation components in Eqs. @) and 



5a = ia5uo d a e * + (iwf 4 - 7 * ot ) 5a + e* b5a) 



+e*a*5b + cfbSJ + ^2^5A in + y^I^Sv™ 



5b 



^/¥c5v s a c + ^f5vt\ 
ibdu$ et + {iuit et - jt ot ) 5b - eaSa 
-\a 2 5e + ^2^5B m + ^2^5v, 



(10) 



out £„,out 
b 



(ii) 



where the intra-cavity and input field vectors are defined 
by 



Xc = 



( 5a \ 

5a^ 
5b 

V W 



Xi n = 



( SAi; 

5Bi, 
V SBl J 



(17) 



the varuum field vectors associated with the output cou- 
pling, absorption loss, and scattering loss are respectively 
defined by 



Vo 



( K nt \ 

Sv outf 
Sv out 

\ K utt 
V H sct J 



\ 


( Svf s \ 








6vf s 


) 


V K 6st / 



(18) 



and the field vector due to photothermal fluctuations, 
X pt , is split into the fluctuating nonlinear coupling con- 
stant component and the fluctuating cavity detuning 
component, 



The coherent components of the equations of motion 
are obtained similarly by taking the coherent terms in 
Eqs. (0J and JSJ) in the steady state, assuming that the 
pump field is undepleted (ea 2 <C \r^B in ), 



= 
= 



2~fl n Bi, 



(12) 
(13) 



where 



X f = 



X p t — X t + X u 



/ 


a*b5e^ \ 


( 




ab*5e 






~\a 2 5e 


\ 




V 



ia5w^ et 
-ia*5w iet ' < 
ib5wf et 



-ib*5w 



det\ 



(19) 



(20) 



from which we find the coherent intra-cavity field ampli- 
tudes, 



^2qfA in [(iu det + 7* ot ) + e* |6| e''^*- 2 *")] 



b si 



ll ot 



la 



B, 



-.detl 



\ef\b\ 2 



ILJ 



det 1 ''' 171 ' 



,(14) 
(15) 



where <j) a and (f>b are the phases of the fundamental and 
second-harmonic input fields such that Ai n = \Ai n \e 1 '^' 1 
and B in = \Bi n \e l ^ b . The relative phase of the funda- 
mental and second-harmonic fields determines whether 
the fundamental field is parametrically amplified or de- 
amplified. 

Eqs. I)10[l. (|ll|l . and their correlated fluctuation oper- 
ators can be rewritten in a compact form, 



- M 0Ut V 



-M SC V SC + M Qbs V abs + X, 



pt ■ 



(16) 



5e, 5uj det , and their adjoints will be derived in the fol- 
lowing section. The coupling matrices associated with 
the intra-cavity field, input coupling, output coupling, 
absorption, and scattering are respectively defined by 

Mr = 





• -det _ tot 
luJ a la 


e*b 


e*a* 









eb* 


^et _ T tot 





ea 






—ea 












V 





— e*a* 





-iuof 1 - 


ll ot 



M in = 



M n 



( V2tF 



V o 










V o 





























out 









VW* 7 



3 



5s 













V 

o 



o 









\ 



V 










o 

y^W J 



In terms of frequency components, defined by |]J 

/oo 
Q(ty nt dt, 
-oo 

Eq. 1|16|) becomes 

iQ# c = M c X c + M in X in + M out V out 



+M SC V SC + M abs V abs + X, 



pi ■ 



(21) 



(22) 



(23) 



where fl is the sideband frequency relative to u a . The 
last field vector X pt will be derived in the next section. 
The commutation relations in Eq. imply 

[s(O), s(O')] = 0, [s(tt),tf(n')] = 6(9, -ft), (24) 



for 



= S A m , SB m , 8v° a ut , 6v% ut , <5< c , <5^ c , 5v% 



b.i 



and 



Sv b bs , and the commutation relations between any two 
different states is zero. 



III. PHOTOTHERMAL NOISE 

Throughout this paper, we consider only type I phase- 
matching which is simple and has been shown to gen- 
erate squeezing at low frequencies. The photothermal 
noise is described as follows: squeezing is degraded by 
fluctuations in (1) the nonlinear coupling strength and 
(2) the cavity resonance frequencies, caused by tempera- 
ture fluctuations due to fluctuations in the photon power 
absorbed in the crystal. Since the nonlinear coupling 
strength is a function of the refractive index along the 
ordinary and extraordinary axes and the crystal length, 
fluctuations in the temperature of the crystal cause fluc- 
tuations in the nonlinear coupling strength. In addition, 
since the cavity resonance frequencies are functions of 
the crystal length at both the fundamental and second- 
harmonic frequencies, fluctuations in the crystal's tem- 
perature cause fluctuations in the resonance frequencies. 
In general, nonlinear crystals are absorptive and there- 
fore, the fluctuations in the crystal's temperature are di- 
rectly coupled with the fluctuations in the amplitudes 
of the input fields. Here we do not consider the three- 
dimensional expansion of the crystal. 



A. Power Absorption in the Crystal 

The fluctuations in the temperature of the crystal are 
due to the fluctuations in the optical power absorbed 



into the crystal, which is directly related to the intra- 
cavity field fluctuations. Taking into consideration that 
the absorption occurs over the entire length of the crystal, 
the total absorbed power is given by 



P, 



b.i 



SP a 



b.i 



abs 



abs 



flLJbBlu.Bn 



abs 



— 1iUa 



\A abs \ 2 + \A abs \(SAl bs + 5A 

abs J 



\B a 



b.i 



\B abs \(6Bl s +SB abs ) , (25) 



where the fluctuation terms are 



SB abs = ^2^ s Sb-Sv b abs , 



(26) 



SB, 



bs 



feifstf-Sv?*, (27) 



and the coherent terms are 



A abs = ^/2 7 f s a, 
B abs = \f^f s b, 



A* — 

A abs — 



b: 



bs 



2^ bs i: 



(28) 
(29) 



Hence, we find the power fluctuation term SP a bs 

SPabs 

= fiLu a \A abs \{SA abs +6Al b3 ) 
+hu b \B abs \(5B abs + 5Bl bs ) 



= fw a ^2^\a\ yf^i&a + 5<f) - (Sv a a bs + Sv^) 



+hu b ^f*\b\ yfi^iSb + Stf) - (Svf s + <5< fcst ) 

(30) 

Assuming that the power absorption in the crystal is 
uniform over the crystal length, it is directly coupled with 
chan ge in the crystal's temperature through the equa- 
tion I2ll 



CpV \5T+^j = SP abs , 



(31) 



where p is the crystal density, V is the mode volume, C 
is the specific heat, and tt is the thermal relaxation time 
of the crystal, tt sets the critical frequency (adiabatic 
limit) for the response to the fluctuations in the optical 
power and is therefore given by \PA 122^ 



T T = 



1 



Cprl 



(32) 



where k is the thermal conductivity of the crystal and 
tq is the radius of the nonlinear interaction between the 
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seed and pump fields, assuming they have a Gaussian 
transverse profile and the interaction distance is within 
the Rayleigh range of the fields. Here, we have assumed 
that the radius of the beams is much smaller than the 
length of the crystal and the cross section of the beams 
is much smaller than the cross section of the crystal. We 
then find the associated temperature fluctuations ST in 
the frequency domain, 



ST 



SPab 



(in + n T )c P v 



(33) 



B. 



Fluctuations in the Nonlinear Coupling 
Strength 



In this section, the fluctuations in the nonlinear cou- 
pling strength are calculated for given temperature fluc- 
tuations. This result will then be used in Eq. 1)23(1 . 

The nonlinear coupling constant e is a function of the 
phase mismatch parameter Ak defined by Ak = 2k a — k b , 



■am Akz 
c = n ze 2 sine , 



(34) 



where kq is a constant. The refractive index of a non- 
linear crystal such as magnesium-oxide doped lithium- 
niobate MgO:LiNb0 3 is dependent on the temperature 
through the photo-refractive effect, which is used for 
achieving type I phase-matching. The temperature and 
wavelength dependence of the phase-matching condition 
for MgO:LiNb0 3 is described by the Sellmeier equa- 
tion [23, which can be approximated around the opti- 
mum temperature Tq at the fundamental frequency 



where £ is a constant whose value depends on the crystal's 
properties, and T is the crystal's temperature. 

The fluctuations in the crystal's temperature cause the 
fluctuations in the nonlinear coupling strength 5c through 
the photo-refractive effect and thermal expansion, 



5c 



dc 



SAk 



dc 



Sz 



dAk dz 
dc dAk dc dz 
dAk~dT~ + ITzdr 
dc \ ^ ( dc 



ST 



dAk 



dz 



ST, 



(36) 



Here we have used the Selmeier equation and the thermal 
expansion equation, 



dAk 

~dT 

dz 

dT 



(37) 

az, (38) 

where a is the linear thermal expansion coefficient. From 
Eq. (EMI) we obtain 



dc 
'dAk 

dc 
dz 



iz 1 z Akz 
" Y~Afc + 2 COt — 



z ( Akz 



(39) 
(40) 



Substituting Eq. into Eq. (JHEJ, we express the 

effect of the fluctuating nonlinear coupling constant in 
terms of X c and V ab s , 



X t = M c £ X c + Mf s V abs , 



(41) 



Ak = aT-T ), 



(35) 



where 



1 

iQ, + tt T 



1 



ab*C a y/2^ db*C a ^f 



db*C bx /2^ bs 



1 =*2 



1^2 



1 = *2 



a 2 C a 
CI 



1 

2 a °Q 



( -a*bC* -d*_bC* 
-ab*C a -db*C a 



d*bCt 



db*C bx /2^ bs 



-a*bCl -a*bC; \ 
-ab*C b -ab*C b 

2° °b 2° °b / 



\d 2 C } 



¥* 2 



(42) 



(43) 



and 

C a 



%( J j a ^/2 1 f s \d\ 



C P V 



TlLJl 



CpV 



( de 



2^ bs \b\ r 



\dAk 

(— 

\dAk 



az 



,(44) 



(45) 



C. Fluctuations in the Cavity Detunings 

In this section, the fluctuations in the optical path 
length are calculated for given temperature fluctuations. 
This result will then be used in Eq. I|23|l . 



The photothermal fluctuations couple to the fluctua- 
tions in the optical path length from the following two 
mechanisms: the photo-refractive effect and thermal ex- 
pansion. The fluctuations in the optical path length 
can be converted into cavity resonance frequency fluc- 
tuations, using [24| 

S< = - 2 -p(^ + a a )sT, (46) 

^ = - 2 ^(l^ + a b )sT, (47) 
\ b \n b dT J 



and therefore, assuming that the laser frequencies are 
stable, the cavity detuning fluctuations are 8ui^ et = 5uj^ 
and 5uj^ et = 5u> b . 

Substituting Eq. (|33|l into Eq. (|46|l . we similarly write 
the effect of the fluctuating cavity detuning component 
in terms of X c and V a bs, 

X e = M^X c + M a J >s Vabs, (48) 

where 



bs _ 



in + n T 



1 

in + n T 



iaK a U ay /2^ idKJ\ a ^2^ ia^n k ^f idK a n b ^2^ b 



i*K a U a ^2X^ ~ia*K a U a ^2^f s -ia*K a U h ^f s -ia*K a U h ^f s 



bs 



ibK b Il a y/2^ 



hs 



bs 



\ - fb*K b u a ^2^ -ib*K b n a ^/¥rs 

I -idKJl a -idKJl a -idKJl b -idK a Il b 

id*_K a Il a id*_K a Il a ia*K a Tl b id*K a Tl b 

-ibK b Il a -ibK b Il a -ibK b Ii b -ibK b Ii b 

\ fb*K b U a fb*K b U a fb*K b U b ib*K b U b 



ibK b n b ^j2^ bs 

-ib*K b tt b 



27i 



abs 



ibK b n b ^j2^ bs 
-ib*K b n b 



2l abs J 



(49) 



(50) 



where 



Ha = 

n 6 = 

K a = 
K b = 



hu a ^2 1 f s \d\ 
C P V 



CpV ' 
27rc / 1 dn a 
A a \n a dT 
2?tc ( 1 dn b 
As \n b dT 



a b 



(51) 

(52) 
(53) 
(54) 



IV. QUADRATURE FIELD AMPLITUDES 
WITH THE PHOTOTHERMAL NOISE 

Now that we have described the equations of motion, 
the absorbed power fluctuations, and the associated fluc- 
tuations in the nonlinear coupling strength and cavity 
detunings, we are in a position to put these equations 
together and solve the equations of motion with the pho- 
tothermal effect. We also discuss a limiting case in which 
the quadrature variances can be approximated to simple 
analytic forms under realistic assumptions. 

Substituting Eqs. (JITJ) and igSJ into Eq. (JTHJ) yields 



and substituting this into Eq. 1)23(1 gives 
(»ni - Mc - Mf - M C JX C 

+ M SC V S( 

+ (M tt6s + Mf s +Mf)V ais , 



(56) 



where I is the identity matrix. We thus find the intra- 
cavity field fluctuations 



X c = (iQl - M c - M° — M£ 



M in X in + M 0Ut V 



-M sc V sc + (M abs +Mf s + M a J> s )V abs . (57) 



Defining the extra-cavity field vector by pjSf 



X n 



( SA out \ 
SB out 



(58) 



X pt = (M c e + M.l)X c + (Mf + MT)Va bs , (55) we find 
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= M out (iOI - M c - M° - M° 
+M out (iOI-M c -M^-M^ 



^Min-Yjn + [M out (im - M c 
- 1 M, c V, c + M out (im-M e - 
I 



-M^-M^IVUt-l] % 



a bs 



MZ bs )V abs .(59) 



It is important to note that as the photothermal effect 
is turned off by setting VtF^ — * an( l \J lb bs ~ ¥ 0, 
the photothermal coupling matrices M e c , M e abs , and 
M* as well as M a b s all become zero, and then Eq. i|58|) 
reduces to the solutions of the field evolution equations 
without the photothermal effect, 



Xout = M out (ifil - M c ) 'MinAi, 



M out (ffiI-M c ) 'Mo 



-Mout (ini - M c ) M SC V SC . 



V out (60) 
(61) 



We define the amplitude and phase quadrature field 
fluctuation amplitudes in the frequency domain relative 
to the fundamental frequency respectively 



8X1 (SI) 



Ss(u a + n) + s§i (<jj a -n), (62) 

i (6S(uj a + ft) - 63* {uj a - fi)) , (63) 



for s = A m , A out , B m ,B out , v° ut , v° b ut , v s a c ,v s b c , vf s , and 
v b bs . The commutation relations (|24|l imply the follow- 
ing values for the commutators of the quadrature field 
amplitudes and their adjoints: 



5Xl,5X 2 s 



5X 2 ,5X 



it 



-2i5(VL - n') (64) 



for s = A m , A out , B in , £ out , v° ut , u 6 out , < c , «g c > < 6s , and 
v b bs , and all others vanish. 

It is convenient to express Eq. I|17|) in terms of the 
quadrature field amplitudes, 



y* 

^out 



AX„. 



out •> £X ~in. ■ L *- cx "tfm 

v; c = AV SC , V* abs = AV absi 



V* t 

Y out 



AV 



(65) 



where 



A = 



/ 1 1 
i -i 
1 f 

V i -i 



and 





/ 5X A out ^ 




f sx A in \ 


y* — 

^out — 


( \ 


x* — 

) -in — 


* X l n 

6X kn 

\ sx kJ 

Ux^\ 


^out 


5X^ t 

K ^ J 

( sZl*. \ 


. K = 


5X% 
5X^ 

V 6x k J 


v\ = 

v abs — 


5X 2 abs 

{ ) 







(66) 



(67) 



Eq. can be rewritten as 

X out = ®inX* n + QoutVout + @s C V* c + Q a b s V* bs , (68) 



where the quadrature field coupling matrices are defined 
by 



Q in = AM mi (fflI-M c -M^M^)' 1 M m A- 1 , 

O out = A[M out (ini-M c -M^-M^,)- 1 M out -l] A" 1 , 

6 SC = AM OTt (iOI - M c - M £ c - M^-'M^A" 1 , 

Q abs = AM^tCim-Mc-M^-M^-^M^ + M^ + Mf^A- 1 . 

I 

Normalized amplitude and phase quadrature variances are given by |19| | 

v s 1 (ri) = (\5xt(n)f) 1 v?w = (\sx:ii>e 



(69) 



(70) 
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for s = A out , B out , A m , B m ,v° a ut ,v° b ut ,v s a c ,v s b c ,vf s , and 
v b bs respectively. The normalized quadrature variances 
of the fundamental output field can be written as a linear 
combination of , v££ , V^l , , , V v \ 2 , V$. , 

~ . a- b a O a 

and V a bB |23 . Since the vacuum fields that couple in at 
the optical losses are in the minimum uncertainty state, 

V 4 ' 2 , = V 1 ' 2 , = Vho = V 1 : 2 = V 1 ' 2 = V 1 '? = 1 (711 

Therefore, we find the normalized amplitude and phase 
quadrature variances of the fundamental output field re- 
spectively, 

= |e£ x) (n)| 2 t>i in (n) + \eg 2) (n)\ 2 vl n (si) 
+|e^ 3) (n)| 2 Fi in (n) + \e^\Q)\ 2 v 2 in (n) 

+£ (|e^(n)| a + |eW)(Q)| 2 + \e ( 2 ] m 2 ) ,(72) 
3=1 

+ £ (|e^(n)| 2 + |e^)(f7)| 2 + \q%J m 2 ) ,(73) 
3=1 

where the superscripts (ij) of 0's denote the matrix ele- 
ments, is the amplitude noise coupling con- 
stant of the seed/pump field, 9^ 2 is the phase 
noise coupling constant of the seed/pump field, and the 



For the amplitude quadrature variance, 



rest of the O's are the amplitude and phase noise cou- 
pling constants of the vacuum fields at the fundamental 
and second-harmonic frequencies. Note that the normal- 
ized quadrature variances are completely characterized 
by the normalized quadrature variances of the two input 
and vacuum fields and the coupling constants. 

If the seed and pump fields are shot-noise limited, 
V A ] = Vg 2 = 1, and the quadrature variances reduce 
to 1 " 

3=1 

+ |eM(fl)| 2 + |e<£ ) (n)| a ), (74) 
3=1 

+ \ei 2 c ^m 2 + \e { a l\n)\ 2 ). (75) 

We now turn to the discussion of a limiting case to 
obtain simple analytic forms by making the following 
assumptions that are applicable to most practical cases 
of squeezing: (1) the pump noise is comparable to the 
seed noise, (2) the fundamental intra-cavity field is much 
weaker than the second- harmonic intra-cavity field (\a\ <C 
|6|), (3) there is no intra-cavity scattering, (4) there are 
no average cavity detunings at both the fundamental and 
second-harmonic frequencies, (5) the frequency of inter- 
est is within the linewidth of the OPA cavity, and (6) we 
are only interested in maximum phase squeezing (</>;, = 0) 
or maximum amplitude squeezing (4n, = tt) so that any 
term proportional to b — b* is zero. Under these assump- 
tions, we can make the following approximations: 



(ii) y/ltr-yrWa * + eb* +?b) 



0(13) 



0(11) 



Q (13) 



e (ii) 

^abs 



o(13) 
^abs 

0(H) 



li° t2 -\e\ 2 \b\ 2 







( 12) u/^^(e*b-eb*) 



li° t2 -\e\ 2 \b\ 2 



1 2V2t^ra Ut lf s ir + ?b) + a*bQ( 7 * ot + eb*)] 

n-ifi T 7 * ot ( 7 * ot2 - |e| 2 |5| 2 ) ' 

- 7 * ot2 + gTTST + l°a ut (fb* + e*b) + |e| 2 |fr| 2 ( i 2) ^ ir a ut (e*b-eb*) 



0(14) „ 



ii° t2 -w\b\ 2 



la 



ton _ |g|2| 6 |2 



1 2V2iJr a ut Jb lbs "f°b Ut + 1*1) + a*6C 6 *( 7 *°* + eb*)] 



Q — iVLt 

A / 7 g«»» 7 g"t(2 7 «° t + a* + rfe) 

li ot2 -\t\ 2 \b\ 2 



il ot {ii° t2 -ww) 



0(14) _ o 



( i 2) ly/jfs^^b-eb*) 







abs 



toll _ |g|2|5|2 



1 i v^F(2 7 f * - 7 < ot ) [ab*C b ( 7 *°* + e*6) + a*bC* b ( 7 * ot + a*)] 



f2 — iVLt 



ii°\ii ot2 - ww) 



0(14) „ o 

O, 



0i 12 » * 0, 



e^) * o, 



(76) 
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For the phase quadrature variance, 



o(21) 



e 



(23) 



e 



@(21) 







(23) 



(21) 



(=)(23) 

^abs 

Q (21) 



7 tot2 _ | e -|2|£|2 > 



e 



/a 



tot2 _ |g|2| fo |2 



1 



Ha ut (e*b-eb*) 



7 *o t(7 t ot 2_| e -|2| 6 | 2) 

,(22) „ - 7 * ot2 + gTTjTT ~ 7a °" t (a* + gj) + \e\ 2 \b\ 2 

I a 



e C") ~ 



«(22) 
^out 



v tot2 „ |g|2|5|2 



1 



2V2 v / 7Q oiI *7f s 7r t [^W* ~ rl) - (flCfa™ - e&*)] 



|2|5|2) 



7 to t 2_| e -|2| 6 |2 



V 

/ a 



to*2_ |e|2|5|2 



1 v/27<p(27 b afcs - 7 * ot ) [a&*C b ( 7 r - rl) - a*~bC* b ( 7 f* - eb*)] 



n — iftT 



7r(7 * ot2 _| e -|2| 6 | 2) 

e(f) ~ o, e(f) ~ o, 



"out — U ! 



(24) _ o 



ei 24 ) ~ o. 



(77) 



Note that in this limiting case, the noise coupling con- 
stants, 613, 0^, 9^, ©™, ®llt, and 9f bs , have the fre- 
quency dependence of 1/(17 — iflr), and therefore they 
increase as the frequency decreases for f2 3> f2x, degrad- 
ing the squeezing level at low frequencies. In Section Ivl 
we will discuss which couping constants are dominant at 
high frequencies, low frequencies, and intermediate fre- 
quencies. 

V. RESULTS 

In this section, we discuss various cases of the in- 
fluence of the photothermal noise on squeezed quadra- 
ture variances in both amplitude and phase quadra- 
tures. The most significant photothermal effects are seen 
in the phase quadrature, and therefore we mainly dis- 
cuss quadrature variances in the phase quadrature. Sec- 
tion IV Al presents such results. Section IV Bl discusses 
the effect of squeezing with the photothermal noise on 
a conventional gravitational wave interferometer when a 
photothermal-noise- limited squeezed field is injected into 
it. The following plots are obtained from the exact nor- 
malized quadrature variances in Eqs. C2J and G3 with 
realistic values for OPA parameters, which are listed in 
Table |I] The effect of green- induced infrared absorption 
(GPJIRA) [28| is not considered in this paper. 

A. Normalized Quadrature Variances with the 
Photothermal Effect 

The amplitude quadrature is relatively immune to the 
photothermal noise for the following reasons. The intra- 
cavity fundamental field is deamplified in the degener- 
ate parametric oscillation, and thus the noise coupling is 



smaller than in the phase squeezing case. Moreover, in 
the ideal case, the system is held on resonance and op- 
erated at the phase-matched temperature. The detuning 
fluctuations do not couple into the amplitude quadra- 
ture for a cavity on resonance since the frequency deriva- 
tive of the amplitude response of the cavity is zero. Fig- 
ure |21 compares the effect of the photothermal noise be- 
tween the amplitude and phase squeezing cases. (Note 
that they are not obtained simultaneously from the OPA; 
different pump phases are required.) The photothermal 
noise is therefore significant in the phase quadrature and 
relatively unimportant in the amplitude quadrature, in 
most practical cases. In the absence of the photothermal 
noise, the normalized quadrature variance would be flat 
within the OPA linewidth. 

As can be seen in Figures |2 El GI and in which 
the normalized quadrature variances versus frequency are 
plotted, the squeezing level is cut off at frequencies be- 
low 10 kHz (depending on the OPA cavity parameters) 
due to the photothermal noise which has 1 / (r2 2 + fi?n) 
roll-off in variance. Parameter values used for the fig- 
ures are summarized in Table [I] The high frequency 
cutoff is due to the linewidth of the OPA cavity be- 
low which the seed field is squeezed. At frequencies 
above the high cutoff frequency, o ^V and Q^ut dom- 
inate in the amplitude and phase quadratures respec- 
tively. At frequencies between the two cutoff frequencies, 

e^.eS.V.e^.e^.e^, and 0<£> dominate (de- 
pending on the OPA cavity parameters) in the amplitude 
and phase quadratures respectively. The photothermal 
cutoff frequency is greater than the adiabatic limit f2^ 
in most practical cases, and therefore, at low frequen- 
cies above Q T , 0<f, 9™, e| 23) , and that have 

the frequency dependence of 1/Q dominate (depending 
on the OPA cavity parameters) in the amplitude and 
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FIG. 2: The comparison of normalized amplitude and phase 
quadrature variances relative to the shot noise vs frequency 
with the photothermal effect. The amplitude quadrature is 
relatively immune to the photothermal noise in most practical 
cases. The seed power is 1 mW. The pump power is 0.5Pth. 
The input fields are shot noise limited. 



FIG. 4: The normalized phase quadrature variance relative 
to the shot noise vs frequency with the photothermal effect 
for different seed powers. The photothermal cutoff frequency 
is higher for a higher seed power. The pump power is 0.5Pth. 
The input fields are shot noise limited. 




FIG. 3: The normalized phase quadrature variance relative 
to the shot noise vs frequency with the photothermal effect for 
different pump powers. The photothermal cutoff frequency is 
higher for a higher level of squeezing whereas it is lower for a 
lower level of squeezing. The seed power is 1 mW. The input 
fields are shot noise limited. 

phase quadratures respectively, as can easily be seen in 
Eqs. (|76|l and (|77|l . At frequencies below the adiabatic 
limit Qt, the quadrature variances become flat. The 
domination of these 0's in each frequency band is valid 
regardless of which quadrature variance is squeezed/anti- 
squeezed. 

Figure |21 shows the normalized quadrature variance of 
a phase squeezed state as a function of frequency for var- 
ious pump powers. As the pump power approaches the 
OPA threshold, the squeezing level at high frequencies 



FIG. 5: The normalized phase quadrature variance relative 
to the shot noise vs frequency with the photothermal effect 
for different pump noise levels. The photothermal cutoff fre- 
quency is higher for a higher pump noise level. The seed 
power is 20 mW. The pump power is 0.5P th- 



increases, but higher pump power also increases the pho- 
tothermal noise contribution. The photothermal noise 
is largest at low frequencies and limits squeezing to oc- 
cur only at higher frequencies where the photothermal 
noise is small. The increase in the photothermal noise as 
the pump power approaches the OPA threshold is also 
attributable to the increase in the fundamental field am- 
plitude via high parametric gain which increases the pho- 
tothermal noise coupling. Squeezing at lower frequencies 
can be acquired at the expense of the broadband level of 
squeezing. 
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FIG. 6: Left: The normalized phase quadrature variance relative to the shot noise level vs pump power with the photothermal 
effect for different frequencies. The seed power is 1 mW. Right: The normalized phase quadrature variance relative to the shot 
noise level vs seed power with the photothermal effect for different frequencies. The pump power is 0.5P th- The input fields 
are shot noise limited in both graphs. 



Figure 0] shows the normalized quadrature variance of 
a phase squeezed state as a function of frequency for var- 
ious seed powers. Since the coupling of the photother- 
mal effect to the quadrature variances is proportional to 
the seed power, the photothermal noise contribution lim- 
its squeezing to higher frequencies as the seed power in- 
creases. Hence, the seed power should be set as small as 
possible to avoid the photothermal noise. However, the 
reduction of the seed power leads to difficulties in ob- 
taining an optical signal for controlling the phase of the 
seed. Therefore, in practice, the seed power should be 
properly chosen such that it optimizes control stability 
and the frequency of interest is above the photothermal 
cutoff frerequency. A crystal with a smaller absorption 
rate can also reduce the pump noise coupling and there- 
fore the photothermal noise. Its effect appears similar to 
the effect of a lower seed power as shown in Figure 0J 

Figure [S] shows the normalized quadrature variance of 
a phase squeezed state as a function of frequency for var- 
ious pump amplitude noise levels. As the pump noise in- 
creases, the overall squeezing level also decreases due to 
the direct coupling of the pump noise to the quadrature 
variances 0> @> @ • At the same time, the photothermal 
noise induced by the pump noise also increases, driving 
up the photothermal noise limited frequency. 

Figure [B] shows the quadrature variance of a phase 
squeezed state as a function of pump/seed power for 
various frequencies. Without the photothermal effect, 
the maximum squeezing would be achieved at the OPA 
threshold. In the presence of the photothermal noise, 
the squeezing level starts to degrade as the pump power 
approaches the OPA threshold. The photothermal noise 
can also be minimized by reducing the seed power. 



In summary, the undesirable consequences of the pho- 
tothermal noise can be minimized by satisfiying the fol- 
lowing conditions: (1) a low seed power, (2) a quiet pump 
field, (3) a crystal with a low absorption rate, and (4) 
squeeze the amplitude quadrature variance rather than 
the phase quadrature variance. 



B. The Effect of the Photothermal Noise on 
Gravitational Wave Interferometers 

One primary purpose of low frequency squeezing is to 
improve the sensitivity of GW interferometers in the GW 
band which is typically 10 - 10,000 Hz 4]. To implement 
it, a low-frequency squeezed field needs to be prepared 
for injection to the dark port of the GW interferometers. 
However, as discussed in Section [V Al the squeezing level 
of phase-squeezed light is limited by the photothermal 
noise at low frequencies, and hence it places an impor- 
tant limit on the use of squeezed light in the GW inter- 
ferometers. 

For a conventional GW interferometer with arm 
lengths L and mirror masses m, the spectral density of 
the GW noise when a realistic squeezed field is injected 
to the dark port is given by 

5(0) = 

^ Q + V ) fa C ° s2{d + $) + ^ Sin2(0 + $) ) ' (78) 
where 

- vSiS (79) 
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TABLE I: OPA Cavity Parameters 



Parameter 


Symbol 


Value 


Units 


Fundamental Wavelength 


A a 


1064 


nm 


Second-Harmonic Wavelength 


A b 


532 


nm 


Reflectivity of Input Coupler at Fundamental Frequency 


R™ 


99.96 


% 


Reflectivity of Output Coupler at Fundamental Frequency 


R° u 


95.6 


% 


Reflectivity of Input Coupler at Second-Harmonic Frequency 


R l b n 


4.0 


% 


Reflectivity of Output Coupler at Second-Harmonic Frequency 


Rl u 


99.96 


% 


Absorption Rate at Fundamental Frequency 




0.1 


%/cm 


Scattering Rate at Fundamental Frequency 




0.02 


%/cm 


Absorption Rate at Second-Harmonic Frequency 


<Jb b3 


4.0 


%/cm 


Scattering Rate at Second-Harmonic Frequency 




0.5 


%/cm 


Crystal Length 


z 


7.5 


mm 


Nonlinear Coupling Strength 


KO 


800,000 


1/m/s 


Phase-Matched Refractive Index 


n 


2.233 


- 


Specific Heat of Crystal 


C 


633 


J/kg/K 


Density of Crystal 


P 


4.648 


g/cm 


Thermal Conductivity of Crystal 


K 


4 


W/K/m 


Radius of Nonlinear Interaction 


ro 


36 


/an 


Phase Mismatch Constant 




749 


1/m/K 


Thermal Expansion Constant in Ordinary Axis 




5 x 10~ 6 


1/K 


Thermal Expansion Constant in Extraordinary Axis 


otb 


5 x 10" 6 


1/K 


Photo-refractive Constant in Ordinary Axis 


dn a /dT 


3.3 x 10~ 6 


1/K 


Photo-refractive Constant in Extraordinary Axis 


dn b /dT 


37.0 x 10~ 6 


1/K 


Temperature Offset 


AT 


0.001 


K 


Cavity Detuning at Fundamental Frequency 


det 





Hz 


Cavity Detuning at Second-Harmonic Frequency 


det 





Hz 



is the noise spectral density of the dimensionless GW 
strain at the standard quantum limit (SQL) for a GW 
interferometer with uncorrelated radiation pressure noise 
and shot noise, 



2{h/I S Q L )l i 



(80) 



is the effective coupling constant that relates the output 
signal to the motion of the GW interferometer mirrors, 
V^(f2),V^(fi) are the amplitude and phase quadrature 
variances of the input squeezed field with a squeeze angle 
9(H) respectively, and 



cot 



= cot 



4 #) 

" 2(W/s Q l)7 4 
_ fF( 7 2 + ft 2 ) 



(81) 



Here 7 is the linewidth of the arm cavities, Iq is the opti- 
cal power to the beam-splitter of the GW interferometer, 
and I sql is the optical power to reach the SQL. 

Figure [3 shows the noise spectral density for the con- 
ventional GW interferometer in various cases when a 
photothermal-noise-limited amplitude-squeezed field is 
injected into the interferometer. The effect of squeez- 
ing with the photothermal noise on the sensitivity of the 



GW interferometer is plotted in the figure. Since the 
amplitude quadrature is relatively insensitive to the pho- 
tothermal noise, we choose the squeeze angle 8 such that 
the amplitude quadrature variance is squeezed and the 
phase quadrature variance is anti-squeezed. The seed 
and pump powers to the OPA cavity are 10 mW and 
0.5P t h- If such a squeezed light field is used for injection 
without a rotation of the squeeze angle, it degrades the 
spectral noise density at about 200 Hz compared with 
the unsqueezed case although it reduces shot noise at 
high frequencies. If a set of two filter cavities is used 
to give the frequency dependent squeeze angle 9(H) such 
that 9(H) = —$(17) 0, the sensitivity is improved at all 
frequencies. This is because the second term in Eq. (f7S|) 
becomes zero and the phase quadrature variance does not 
couple into the sensitivity curve. However, it is likely to 
be difficult and costly to implement such a squeeze angle 
rotation since the length of the filter cavities is required 
to be long (on the order of kilometers) in order to mini- 
mize losses that destroy squeezing in the process 0, |2t| . 
If a squeeze amplitude filter is instead used before inject- 
ing the squeezed light into the GW interferometer [30| 
such that for a filter linewidth jf, 



S(H) 



h 2 

"■SQL 

277 



Ci(^i+^ 2 V r 2 ) + C2(l + ?7 2 ) 



(82) 
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Frequency (Hz) 



FIG. 7: The spectral noise density, normalized by the 
standard quantum limit (SQL), for a conventional GW 
interferometer with (i) no squeezed input (Unsqueezed), 
(ii) squeezed light injected without the photothermal noise 
(Frequency-Independent), (iii) squeezed light injected with 
a frequency-dependent squeeze angle and the photothermal 
noise (Frequency-Dependent + Photothermal), (iv) squeezed 
light injected with the photothermal noise (Frequency- 
Independent + Photothermal), and (v) amplitude- filtered 
squeezed light injected with the photothermal noise (Ampli- 
tude Filter + Photothermal). The filter linewidth is 2n x 400 
Hz. The input squeezed source is chosen to be photothermal- 
noise-limited amplitude squeezed light since the amplitude 
quadrature is relatively immune to the photothermal noise. 
The seed and pump powers are 10 mW and 0.5Pth respec- 
tively. The seed and pump fields are shot noise limited. The 
anti-squeezed phase quadrature variance has the photother- 
mal cutoff frequency at about 1 kHz and the adiabatic limit 
at about 100 Hz. 



where 



Ci(n) 



C 2 (0) 



(83) 



the phase quadrature variance containing the photother- 
mal noise can be reduced at frequencies below 200 Hz 



although the level of squeezing at frequencies above 200 
Hz is slightly decreased. Here 7/ = 2n x 400 Hz is used. 
Note that a photothermal-noise-limited phase-squeezed 
input field with similar experimental parameters (the 
seed power = 1 mW, the pump power = 0.5Pth) does 
not enable quantum noise reduction below 400 Hz, even 
if the optimal frequency dependent squeeze angle rota- 
tion is applied. 



VI. CONCLUSIONS 

We have derived and solved the field evolution equa- 
tions in the degenerate optical parametric amplifier 
(OPA) with the photothermal noise through the photo- 
refractive effect and thermal expansion of nonlinear crys- 
tals. We also have discussed various cases about the ef- 
fect of the photothermal noise on amplitude and phase 
quadrature variances. We have found that the photother- 
mal noise in the OPA introduces a significant amount of 
noise on phase squeezed beams, making them less than 
ideal for low-frequency applications such as GW inter- 
ferometers, whereas amplitude squeezed beams are less 
sensitive to the photothermal noise and may provide a 
better choice for low frequency applications. This prob- 
lem can be solved by reducing the seed power and pump 
noise and using a nonlinear crystal with a low absorption 
rate in order to decrease the photothermal noise. 
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